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ABSTRACT 

High-resolution observations show that oscillations and waves in prominence 
threads are common and that they are attenuated in a few periods. In addition, 
observers have also reported the presence of material flows in such prominence 
fine-structures. Here we investigate the time damping of non-leaky oscillations 
supported by a homogeneous cylindrical prominence thread embedded in an un- 
bounded corona and with a steady mass flow. Thermal conduction and radiative 
losses are taken into account as damping mechanisms, and the effect of these 
non-ideal effects and the steady flow on the attenuation of oscillations is assessed. 
We solve the general dispersion relation for linear, non-adiabatic magnetoacous- 
tic and thermal waves supported by the model, and find that slow and thermal 
modes are efficiently attenuated by non-adiabatic mechanisms. On the contrary, 
fast kink modes are much less affected and their damping times are much larger 
than those observed. The presence of flow has no effect on the damping of slow 
and thermal waves, whereas fast kink waves are more (less) attenuated when they 
propagate parallel (anti-parallel) to the flow direction. Although the presence of 
steady mass flows improves the efficiency of non-adiabatic mechanisms on the at- 
tenuation of transverse, kink oscillations for parallel propagation to the flow, its 
effect is still not enough to obtain damping times compatible with observations. 

Subject headings: Sun: oscillations - Sun: magnetic fields - Sun: corona - Sun: 
prominences 



1. Introduction 



Prominences and filaments are large-scale magnetic structures embedded in the solar 
corona, and whose plasma density and temperature are akin to those of the chromosphere. 
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High-resolution images of solar filaments (e.g., iLin et al.ll2003l . |2005| . 120071 ) clearly show the 
existence of horizontal fine-structures within the filament body. This observational evidence 
suggests that prominences are composed of many field-aligned threads. These threads are 
usually skewed with respect to the filament long axis by an angle of 20° on av erage, although 
their orientation can vary significantly within the same prominence (jLinll2004l ). The observed 
thickness, d, and length, /, of t hreads are typic ally in the ranges 0.2arcsec < d < 0.6arcsec 
and Sarcsec < / < 20arcsec (ILin et al.l l2005l ). Since the observed thickness is close to 
the resolution of present-day te lescopes, it is likely that even thinner threads could exist. 



According to some models (e.g., Ballester fc Priestfll989l ). a thread is believed to be part of 



a larger magnetic coronal flux tube which is anchored in the photosphere, with denser and 
cooler material near its apex, i.e. the observed thread itself. However, the process that leads 
to the formation of such structures is still unknown. 



There are many evidences of small-amplitude waves and oscillatory motions in 



cent prominences (th is topic has been reviewed by lOliver fc Ballesterl l2002l : iBallesten 12006 



quies; 



Banerjee et al.l 120071 ). Focusing on prominence fine-structures, some observers have de 



tected oscillations and trayeling wayes in individual threads or groups of threads (jYi et al. 



199ll : lYi fc Engvoldlll99ll : iLinI 12004 iLin et all 120071 ) . with periods typically between 3 and 
20 min. In addition, mass flows along fil ament threads, wi t h a flow velocity in t he range 
5-25 km s~^, have been also observed (Zirker et al. 



1998 



Lin et al.ll2003l . 120051 ). More- 



over, it is noticeable that IZirker et al.l (119981 ) and iLin et al. J2003h have detected flows in 



opposite directions within adjacent threads, a phenomenon known as counter-streaming. On 
the other hand, some o bservational works have suggested signatures of wave dam ping in 
prominen ce oscillations (jLandrnan et al. 1977 : Tsubaki fc Takeuchi 1986 ; Lin 2004), but to 
date only Molowny-Horas et al. (1999) and Terradas et al. 120021 ) have studied in detail this 
phenomenon, in particular in two-dimensional Doppler velocity time series. This analysis 
showed that oscillations detected in large areas of a quiescent prominence were attenuated 
after few periods. Although this quick attenuation seems to be a common feature of promi- 
nence oscillations, unfortunately no similar observational study focusing on the attenuation 
of individual thread oscillations has been performed yet. 

From the theoretical point of view, the usual interpretation of thread oscillations is 
in terms of the adiabatic magnetohydrodynamic eigenmodes supporte d by the thread body . 
The first investigation of individual thread v ibrations was perf ormed by lJoarder et al.l (119971 ). 
whose work was extended and corrected by iDfaz et al.l (120011 ) , considering a non-isothermal 
Carte sian thread surrounded by the coronal medium (based on the model by lBallester fc Priest 
19891 ). in the (3 = approximation. These authors found that only the low-frequency os- 
cillatory modes are confined within the dense region, and that perturb ations can a chieve 
large amplitudes in the corona at long distances from the thread. Later, iDiaz et al.l (120031 ) 
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assumed the same geometry, but took longitudinal propagation into account, and obtained 
a better confinement for the per t urbat ions. Considering a more realistic and representative 
cylindrical geometry, iDiaz et al.l (120021 ) found that a non-isothermal cylindrical thread sup- 
ports an even smaller number of trapped oscillations and that perturbations are much more 
efficiently confined within the cylinder in comparison with the Cartesian case. It is worth to 
m ention that the co llecti ve oscillations of multithread systems have also been in vestigated 
bvlDiaz et al.l (120051 ) and iDiaz fc Roberts! (120061 ) , again in Cartesian geometry. See iBallester 
(120061 ) for a review about theoretical works. 



The effect of steady mass flows on the oscillatory modes of magnetic structures has been 
the oretically investigated in some works. The most relevant ones for the present investigation 



are lNakariakov &: Roberts! (119951) . who studie d the effect of steady flow in coronal and photo- 
spheric slabs, and !Terra-Homem et al.! (12003! ). who extended the former study to cylindrical 
geometry. In addition to producing a shift of the oscillatory frequency, both papers show 
that the main effect of the flow is to break the symmetry between parallel and anti-parallel 
wave propagation to the flow direction and, for sufficiently strong flows, slow modes can only 
propagate parallel to the flow direction, anti-parallel propagation being forbidden. 

Turning to the damping of oscillations, its theoretical investigation has been under- 
taken by a numb er of recent papers. Among the proposed damping mechanisms to explain 
the attenuation (IBallai! 12003! ). non-adiabatic effects are the most ext en sively investigated 
to da te, although other candidat es like wave leakage (ISchutgens Il997al jbl: I Schuteens k. TothI 
1999! ) and ion-neutral collisions (IForteza et al.ll2007! ) have been also studied. By removing 
the adiabatic assumption and taking into account thermal conduction and radiative losses 
as damping mechan isms, some works hav e studied the time damping in a homogeneous. 



unbo unded plasma (jCarbonell et al.ll2004j ). in an isolated promin ence slab (iTerradas et al. 



20051 ). and in a prominence slab embedded in the solar corona (jSoler et al.l l2007al |2008| ). 
The common main result of these studies is that non-adiabatic mechanisms can explain the 
observed damping times in the case of slow modes, whereas fast modes are much less atten- 
uated by non-adiabatic effects. It is important to note that all these referred articles have 
studied the wave attenuation in models which attempt to represent the whole prominence 
body. Thus, none of them have neither considered the prominence fine-structure nor mass 
flows. 



More recently, ICarbonell et al.l (120081 ) have performed the first attempt to study the 
combined effect of both non-adiabatic mechanisms and steady flows on the time damping of 
slow and thermal waves in a homogeneous, unbounded prominence plasma. These authors 
found that the mass flow does not modify the damping time of both slow and thermal waves 
with respect to the case without flow, but the period of the slow wave increases dramatically 
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for flow velocities close to the non-adiabatic sound speed. Moreover, the thermal disturbance 
behaves as a pr opaKating mode i r i the p resence of flow. The present work goes a step forward 
with respect to ICarbonell et al.l (120081 ) since a more complicated geometry is assumed here. 
Our aim is to describe the effect of both mass flow and non-adiabatic effects on the oscilla- 
tions supported by an individual prominence thread modeled as a homogeneous and infinite 
cylinder embedded in an unbounded and also homogeneous corona. For simplicity, the hot, 
coronal part of the magnetic tube that contains the thread is not taken into account. Gravity 
is also discarded. The inclusion of gravity would involve the consideratio n of a more compli 



cated magnetic structure able to support the prominen ce material (e.g.. iBallester &: Priest 
19891 : ISchmitt &: Degenhardtl Il995l : iRempel et al.lll999l ). However, such kind of configura- 
tions are unstable and incomplete, since they do not incorporate the physics that lead to 
stable solutions and, moreover, the obtained prominence widths are much smaller than those 
observed. Although some prominences and threads show an unstable behavior, ther e is a 
number of observations of stable oscillating threads (e.g., iLin et al.ll2003l . |2005| . 120071 ). For 
this reason and since the present work is focused on the study of waves in stable threads, 
we neglect the effect of gravity and consider a stable simplified model that allows us to deal 
with wave solutions. 

This paper is organized as follows. The description of the model configuration and the 
basic equations for the discussion of linear non-adiabatic waves are given in § |2j Then, the 
results are presented in § [3], first for the case without flow and later extended by including 
a steady mass flow in the equilibrium. Finally, § S] contains our conclusions. 



2. Model equations 

The model configuration considered in the present work (Fig. [1]) is made of a homoge- 
neous and isothermal plasma cylinder of radius a with prominence conditions (density pp 
and temperature Tp) embedded in an unbounded corona (density pc and temperature T^). 
The cylinder is also unlimited in the axial direction. The internal temperature and density, 
as well as the coronal temperature, are considered as free parameters. However, the value 
of the coronal density is given by imposing total pressure continuity across the interface 
between the flux tube and the external medium. In all the following expressions, a subscript 
indicates local values, while subscripts p and c denote quantities explicitly computed with 
prominence and coronal parameters, respectively. 




Fig. 1. — Sketch of the modeL 
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Given the model geometry, we use cylindrical coordinates, namely r, and z for the 
radial, azimuthal and longitudinal coordinates, respectively. The magnetic field is uniform 
and orientated along the cylinder axis, Bq = -Bq^z, -Bo being the same constant in the thread 
and in the coronal medium. A steady mass fiow is assumed along the ^-direction, whose fiow 
velocity can be different in the cylinder and in the corona. Thus, Up = f/pe^ and Uc = t/cCz 
correspond to the steady flow in the flux tube and in the corona, respectively. 



Parallel thermal conduction to the magnetic fleld, radiative losses, and heating are 
considered as non-adiabatic effects. We assume that the plasma is fully ionized, and so the 
cross-fleld or perpendicular thermal conduction is absolutely negligible. The contribution of 
neutrals to the ther mal conduction in a partially ionized plasma has been investigated by 
Forteza et al.l (l2008l . in preparation). Radiation and heating are evaluated together by means 
of the heat-loss function, L{p , T) = x*pT°' — hp°'*T^* , in which radiation is parametrized 



h* teosner et al. 



1978 



wit h Y* and a (Hildnerl Il974f) and the heating sc enario is given by the exponents a* and 



Dahlburg fc Mariskal Il988l ) . Since the result s for different heating 



scenarios do not show signiflcant differences f or prominence conditions (jCarbonell et al.ll2004 



Terradas et al.ll2005l : ISoler et al.ll2007al . l2008l ). here we restrict ourselves to a constant heating 
per unit volume (a* = 6* = 0). 

The basic equations for the discussion of non-adiaba t ic osci llations supported by this 
equilibrium correspond to equations (l)-(6) of ISoler et al.l (l2007al hereafter Paper I), whose 
linearized version applied to the present case is given next. 



dpi ^ dpi 
dt ° dz 



PoV ■ vi = 0, 



9vi dwi 



dt 



dz 



-Vpi + - (V X Bi) X Bo, 
P 



(1) 



(2) 



9pi . jj dpi 2 (dpi . 



PqUq 



dpi 
dz 



+ (7-1 
dBi 



PoLpPi + PqLtTi - K\ 



d^Ti 
dz"^ 



(3) 



dt 



dBi 



V X (vi X Be 



dz 

V-Bi = 0, 

Pi _ Pi _ 7i 
Po Po ^0 

Here po, po, Tq, Bq, and Uq are the equilibrium gas pressure, density, temperature, magnetic 
field, and fiow velocity, respectively. On the other hand, pi, pi, Ti, Bi = Brer + B^e^ + BzGz^ 



0. 



(4) 
(5) 
(6) 
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and vi = Vrer + v^eip + VzGz are the linear gas pressure, density, temperature, magnetic field, 
and velocity perturbations, respectively. Finally, Cg = ^ is the adiabatic sound speed 
squared, = 10~^^T^/^ Wm~^ is the thermal conductivity parallel to the magnetic 
field, and Lp and Lt are the partial derivatives of the heat-loss function with respect to 
density and temperature, respectively (see Paper I for details). 

Assuming perturbations of the form f\{r, if) exp {iut — ik^z), where uj is the (complex) 
oscillatory frequency and the (real) longitudinal wavenumber, one can combine equa- 
tions (IH]) and ([6]) to obtain the following relation between the perturbed pressure and density. 



Pi = Kpu 



(7) 



with 



IVq — 

7 



(7 - 1) ( f + _ ^ \ + 



PO 



PO 



where Qq = u — U^k^ is the Doppler-shifted frequency (ITerra-Homem et al.ll2003l ). and Up 
and ujt are defined in Paper I. We see than the complex quantity Aq is a generalization, 
due to the presence of the fiow, of the non-adiabatic sound speed defined in Paper I. The 
real part of Aq plays the role of the sound speed when non-adiabatic effects are present. By 
means of this definition, one can see that the effect of non-adiabatic terms is to modify the 
medium sound speed, and so they most probably affect slow modes since they are mainly 
governed by acoustic effects. See § 13.2.31 for more details about this quantity. 



Now, following iTerra-Homem et al.l (120031 ). we combine the basic equations and arrive 
to the following expressions. 



0, 



where v\ 



T^-.i|.|r^o 

^ is the Alfven speed squared, T is a linear operator defined as follows, 

dt dz' 



(9) 
(10) 



(11) 



and A and F are the divergence and the z-component of the rotational of the velocity 
perturbation, respectively, 

A = V-vi, (12) 

F = (Vxvi)-e,. (13) 
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Equation ffTOj) is the same as equation (14) of iTerra-Homem et al.l (120031 ) and governs tor- 
sional, Alfven waves, which are not damped by non-adiabatic mechanisms and so are not 
considered in the present investigation. On the other hand, equation repres ents fast an d 
slow magnetosonic waves, together with the thermal or condensation mode (iFieldl Il965l ). 
If non-adia batic terms are neglected, An = Cg and then our equation ([9]) reduces to equa- 
tion (13) of iTerra-Homem et al.l (120031). whi c h in t he absence of flow {Uq = 0) is equivalent 
to the well-known equation (17) of iLighthilll (jl960l ). 



Next, the cylindrical symmetry of the model allows us to write the divergence of the 
velocity perturbation in the following form. 



A = R{r) exp [iut + irup — ik^z) 



(14) 



where n is an integer that plays the role of the azimuthal wavenumber. Expressions for the 
perturbed quantities as a function of A can be found in Appendix \M Now, applying this 
last expression to equation one finds that R{r) satisfies the well-known Bessel equation 
of order n, 

^d^R dR 



dr^ 



+ r— + (mlr^ - n'^) R = 0, 
dr ^ 



with 



■rrin 



7;2 A2 



(15) 



(16) 



(17) 



The radial wavenumber squared, mQ, is in general a complex quantity, hence no pure body- 
like or surface-like waves are possible in non-adiabatic magnetohydrodynamics. If one as- 
sumes that |3fJ(mQ)| > 13^(7710)1 (i-^-' non-adiabatic effects produce a small correction to the 
adiabatic wave modes), the dominant wave character depends on the sign of 3?(mQ). Thus, 
oscillations are mainly body-like if 5J(mo) > and solutions of equation (fT5|) are Bessel func- 
tions. On the contrary, if 3?(mo) < oscillations are mainly surface-like (or evanescent) and 
solutions of equation (fTSll are modified Bessel functions. In this work we assume no wave 
propagation in the coronal medium, so the evanescent condition in the corona is imposed on 
the perturbations, namely 3fJ(m^) < 0. On the other hand, the present ordering of sound and 
Alfven speeds does not permit the existence of surface waves within the fibril, so 3fJ(mp) > 
is assumed. Then R{r) is a piecewise function. 



R{r) 



AiJnijnpr) if r < a 
A2Kn{ncr) if r > a 
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with nl = —ml, Ai and A2 being complex constan ts. J„, and Kn are the usual Bessel and 



modified Bessel functions of order n, respectively (jAbramowitz &: Stegun 1972 ). In order 
to obtain the dispersion relation that governs the behavior of wave modes, we impose the 
continuity of the Lagrangian radial displacement, Vr/fio, and the total pressure perturbation, 
PTi, at the cylinder edge, r = a. After some algebra, the following expression is obtained. 



Ac 



m, 



'Jn (mpa) 



KL 



rira] 



rir 



Kn {ricO) 



(19) 



where ' = d/dr. Equation flTIJl) is formally identical to equation (21) of iTerra-Homem et al. 



( I2OO3I ) since our non-adiabatic terms are enclosed in the definition of nip and Uc- If both non- 
adiabatic effects and the flow are dropped, eq uation (fT9!) simply reduces to the well-known 
dispersion relation of Edwin fc Roberts! (11981: 



The solution of equation (IT^ for a real is a complex frequency, = c<jr + iuji. The 
oscillatory period (P), damping time (td), and the ratio of both quantities are computed as 
follows. 



27r 
I^rI 



= — , 



P 



1 kRi 

27r oji 



(20) 



3. Results 

3.1. Configuration without flow 

In this section, we flrst perform a study of the solutions of the dispersion relation 
(eq. [in]) in the absence of steady flow. In this case, Up = Ur = and so Qp = = 



This situation corresponds to the case studied by lEdwin fc RobertsI (1l983l ) with the addition 
of non-adiabatic effects. Unless otherwise stated, the following equilibrium parameters are 
considered in all computations: Tp = 8000 K, pp = 5 x 10~^^ kg m~^, = 10® K, = 
2.5 X 10~^^ kg m~'^, i?o = 5 G, and a = 30 km. Thus, the characteristic speeds of the internal 
and external media are: Cxp = 11.56 kms~^, Cgp = 11.76 kms~^, fAp = 63.08 kms~^, Ctc = 
163.51kms~^, Cgc = 166.33 kms~^, and v^c = 892.06 kms~^. In addition, both prominence 
and coronal plasmas are taken as optically thin (see Table I in Paper I for the values of 
parameters x* aiid a of the heat-loss function). 

In the absence of flow, the complex oscillatory frequencies obtained by solving equa- 
tion (IT^ for a flxed, real, and positive appear in pairs, Ui = u^ + iui and UJ2 = —u^i + iuji. 
The solution uji corresponds to a wave propagating towards the positive z-direction (parallel 
to magnetic fleld lines), whereas 002 corresponds to a wave that propagates towards the neg- 
ative 2;-direction (anti-parallel to magnetic field lines). For short, we call them parallel and 
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anti-parallel waves, respectively. Both parallel and anti-parallel wave modes are equivalent 
and show exactly the same physical properties when no flow is considered. For the sake of 
simplicity, the results presented in this section correspond to parallel waves, which have a 
positive phase speed. Equivalent results for anti-parallel waves are deduced by considering 
negative phase speeds. 



3.1.1. Dispersion diagram and eigenf unctions of magnetoacoustic modes 



Magnet oacoustic modes supported by a magnetic cylinde r have been extensively inves- 
tigated fe.g.. lSpruitlll982l : Isdwin fc RobertslllQSsI : ICallvUlQSGf ). Fast oscillations with n = 0, 
n = 1, and n >2 correspond to sausage, kink, and fluting or ballooning modes, respectively. 
On the basis that thread oscillations are observed in Doppler time series, this work is mainly 
focused on kink modes, which produce displacements of the thread axis from its original 
position. The fundamental fast kink mode is trapped for realistic values of the thickness and 
width of threads. Regarding slow modes, the fundamental modes and their harmonics are all 
trapped for any value of kz and n, but all of them have an almost identical frequency, as in 
the slab case (Paper I). Thus, we also restrict ourselves to the fundamental slow mode with 
n = 1 for simplicity. An additional solution of equation f|T9|) is the thermal mode which, in 
the absence of flow, has a purely imaginary frequency (see § 13.1.31) . 

Figure [2] displays the phase speed diagram corresponding to the fundamental modes 
with n = and n = 1 (compare this with Fig. 2 of Paper I). One can see that the behavior 
of slow modes and the fast sausage mode is similar to that in a slab. An important difference 
with the results of Paper I is that in the cylindrical case the fast kink mode does not couple 
with the external leaky slow modes enclosed in the region 3?(ctJ < i^r/^z < 5?(Ac), because 

its phase speed in the long-wavelength limit is Ck < 5? (ctJ, with = ^ {p'+p') ^ • 

Next, we plot in Figure [3]the eigenf unctions corresponding to the radial and longitudinal 
velocity perturbations, Vr and Vz, and the total pressure perturbation, Pti, for the funda- 



mental slow and fast kink magnetoacoustic modes. Contrary to the slab case of ISoler et al. 



(l2007bl . Fig. 4), perturbations are efficiently confined within the cylinder for any value of 
the longitudinal wavenumber. Thus, this suggests that the infiuence of the corona on the 
damping of oscillations could be of smaller importance than in a magnetic slab. On the other 
hand, the expected velocity polarization is obtained, the slow mode being mainly polarized 
along the longitudinal direction {vz ^ fr) and the fast kink mode being responsible for 
transverse, radial motions (f^ ^ Vz). 
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Fig. 2. — Phase speed {ooji/kz) versus the dimcnsionless longitudinal wavenumber (fc^a) for 
the fundamental magnetoacoustic modes with n = 1 (solid line) and n = (dashed line). 
The shaded zones are projections of the leaky regions on the plane of this diagram. Note 
the cut-off frequency of the fundamental fast sausage mode. The vertical axis is not drawn 
to scale. 



- 12 - 




Fig. 3. — Modulus of the eigenfunctions (in arbitrary units) corresponding to the radial 
velocity perturbation, Vr-, the longitudinal velocity perturbation, f^, and the total pressure 
perturbation, Ptu as functions of the dimensionless distance to the cylinder axis for the fun- 
damental slow (top panels) and fast kink (bottom panels) modes. The solid line corresponds 
to kza — 1 while the dashed line corresponds to /c^a = 10~^. Note that is not continuous 
at the cylinder edge, whose location is denoted by a vertical dotted hne. 
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3.1.2. Damping times of fast kink and slow waves and comparison with a longitudinal slab 

We now compute the period, damping time, and their ratio for the fundamental fast kink 
and slow modes for a wide range of kz (10~^°m~^ < k^ < 10~^m~^). This range includes 
the values corresponding to the observed wavelengths in prominences (10~^ m^^ ^ ^2 ^ 
10^^ m^^). In Figure m these results are compared with those obtained for a longitudinal 
slab whose half- width is equal to the cylinder radius (i.e., the case analyzed in § 4.4 of 
Paper I). We see that the results of the slow mode are almost identical in both cylindrical 
and slab geometries, hence the reader is referred to Paper I for a description of the slow 
mode behavior. 
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Fig. 4. — Period (left), damping time (center), and ratio of the damping time to the period 
(right) versus the longitudinal wavcnumber for the fundamental oscillatory modes in the 
absence of flow. Upper and lower panels correspond to the slow and fast kink modes, 
respectively. Solid lines are the solutions of the present, cylindrical equilibrium whereas 
dotted lines arc the results of a longitudinal slab (Paper 1). Both curves overlap in the upper 
panels. Shaded zones indicate the range of observed wavelengths. 
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On the contrary, the results of the fast kink mode show significant differences between 
the slab and the cylinder. For small and intermediate kz, both the period and damping time 
in a cylinder are larger than in a slab. This effect is specially noticeable for the damping 
time in the observed range of wavelengths, since its value in the cylindrical case is more than 
three orders of magnitude larger than in the slab. This causes the ratio td/P to be in the 
range 10^ ^ tb/P ^ 10^ for the observed wavelengths. Therefore, non-adiabatic effects are 
much less efficient in damping the fast kink mode in a cylinder than in the slab geometry. In 
addition, we see that the fast kink mode damping time due to non-adiabatic effects is much 
larger than typical lifetimes of filament threads and prominences. Hence, non-adiabatic fast 
kink waves are in practice undamped when these results are applied in the solar context. 
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3.1.3. Thermal mode 



Now we turn our attention to the thermal or condensation mode. The condensation 



i nstabihty has been studied in uii iform, unbounded plasmas ( Field! 1 196 51) . in coronal slabs 
(Ivan der Linden fc Goossenslll99ll ). and in coronal cylinders (jAnlll984l ). Since the thermal 
mode has a purely imaginary frequency, we now assume uj = is, where s is real and often 
called the damping (or growing) rate. The situation s > corresponds to a damped thermal 
mode, whereas s < occurs if the mode is thermally unsta ble. The sign of s can be estimated 
a priori by considering the stability criterion provided by iFieldl (119651 ). 



'^Wp^I + Pp 



T P^T 
J^Tp - TfrJ^pp 



> 0. 



(21) 



Since for prominence conditions this inequality is verified for any real value of kz, the ther- 
mal mode is always a darnped so lution and therefore we expect s > 0. According to 
van der Linden fc Goossend (Il99ll ). see Appendix [B] for more details, the evanescent as- 



sumption in the corona (m^ < 0) together with the body-wave assumption within the fibril 
(rrip > 0) is satisfied in a narrow range of s. However, there is an extremely narrow range 



of kz (1.68 X 10- 



m 



< kz < 1.70 X 10 



-''m ^) in which rrip > and 



> 0, and so the 



evanescent assumption is not verified. Then, the thermal mode does not exist as a non- 
leaky solution in such a range of kz- Despite this, the fundamental thermal mode and all its 
harmonics have an almost identical damping rate s, whose value is also almost independent 
of the azimuthal wavenumber, n. For this reason and for the sake of simplicity, we again 
restrict ourselves to solutions with n = 1 and focus on the fundamental mode. 

Figure [5] displays the spatial distribution of perturbations Vr, Vz, and Pti corresponding 
to the fundamental thermal mode with n = 1. We see that the velocity field is dramati- 
cally polarized along the cylinder axis and that the eigenf unctions are very similar to those 
obtained for the slow mode (compare with Fig. [3], upper panels). On the other hand, the 
damping time (td = l/s) is plotted in Figure [6] as a function of kz. One can see that this 
mode is very quickly attenuated and that radiative losses from the prominence plasma are re- 
sponsible for the attenuation in the observed wavelength range, whereas prominence thermal 
conduction is only relevant for large kz. Coronal mechanisms have a negligible effect. 
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Fig. 5. — Same as Fig. [3] but for the fundamental thermal mode with k^a = 1. 



- 18 - 




10''° 10'' 10'" 10'" lO''^ 

(m'') 



Fig. 6. — Damping time versus the longitudinal wavenumber for the fundamental thermal 
mode. Different linestyles represent the omitted non-adiabatic mechanism: all mechanisms 
considered (solid line), prominence conduction eliminated (dotted line), prominence radiation 
eliminated (dashed line), coronal conduction eliminated (dot-dashed line), and coronal radi- 
ation eliminated (three dot-dashed line). The vertical line with symbols inside the shaded 
zone indicates the range of kz in which the thermal mode does not exist as a non-leaky 
solution (see Appendix [B] for details). 
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3.2. Effect of steady flow 

Hereafter, we include a steady mass flow in the model in order to assess its influence on 
the previously described oscillatory modes. With no loss of generality, we assume no flow in 
the external medium, i.e. Uc = 0. Moreover, the internal flow is assumed flowing towards the 
positive z-direction, i.e. Up > 0. Note that any other possible configuration can be obtained 
by means of a suitable election of the reference frame. 



3.2.1. Phase speed shift 



The reader is referred to iTerra-Homem et al.l (120031 ) for a detailed description of the 
modification of the phase speed diagram due to the presence of flow. In short, the symmetry 
between parallel (wr > 0) and anti-parallel (cur < 0) waves is broken by the flow. Phase 



speeds of slow and fast parallel waves are now in the ranges 



3f^(5Tj+f/p,3f^ Ap +f/, 



and 



[vAp + Up,VAc\i respectively. On the other hand, phase speeds of slow and fast anti-parallel 



waves lie now within the regions 



and [-vxc, -v^p + f^p] 



-3? (^Apj +f/p,-3ft (£Tp)+f/p 
respectively. For a fl ow velocity larger than the internal non-adiabatic sound speed (see 
Carbonell et"aD boosh the phase speed of anti-parallel slow waves is dragged to positive 
values, a nd so they become para ll el wav es in practice. These solutions were called backward 
waves by iNakariakov fc Roberts! (119951 ). An equivalent phenomenon can also occur for fast 



waves for a superalfvenic flow, i.e. f/p > fAp- Note that superalfvenic flows seem to be 
unrealistic in the light of observations. 

Regarding thermal modes, the real part of their frequency now acquires a positive value, 
and their phase speed is equal to the flow velocity. Thus, thermal modes behave as prop- 
agating parallel waves with respect to the static, external reference frame. Although this 
result could be relevant for the observational point of view, since t hermal modes might be 
detected as propagating waves in filament threads (ILin et al.l 120071 ). their extremely quick 
attenuation makes them undetectable in practice. 
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Fig. 7. — Period (left), damping time (center), and ratio of the damping time to the pe- 
riod (right) versus the flow velocity for the fundamental oscillatory modes with k^a = 10~^. 
Upper, mid, and lower panels correspond to the slow, fast kink, and thermal modes, respec- 
tively. Different linestyles represent: parallel waves (solid line), anti-parallel waves (dashed 
line) , and solutions in the absence of flow (dotted hne) . 
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3.2.2. Influence on the damping time 



Regarding the effect of tlie ffow on tlie damping time of oscillations, we plot in Figure [7] 
the dependence of the period, damping time, and their ratio as a function of the flow ve- 
locity. The longitudinal wavenumber has been fixed to /c^a = 10^^, which corresponds to 
a value within the observed range of wavelengths. The flow velocity is considered in the 
range Okms^^ < Up < 30kms~ \ that corresporids to t he observed flow speeds in quiescent 
prominences. In agreement with lCarbonell et al.l (120081 ). the anti-parallel slow wave becomes 



a backward wave for f/p 8.5 km s~^, which corresponds to the non-adiabatic sound speed 
(eq. [S]). This causes the period of this solution to grow dramatically near such flow veloc- 
ity. However, the period of both parallel and anti-parallel fast kink waves is only slightly 
modified with respect to the solution in the absence of flow, and the thermal wave now has 
a finite period, which is comparable to that of the parallel slow mode. 

On the other hand, we see that t he damping time o f both slow and thermal modes is not 
affected by the presence of flow, as in ICarbonell et al. J2008h . Nevertheless, the attenuation 
of the fast kink mode in the present case is influenced by the flow. The larger the flow 
velocity, the more attenuated the parallel fast kink wave, whereas the opposite occurs for 
the anti-parallel solution. This behavior can be understood with the help of Figure [8], which 
displays the phase speed of the parallel fast kink mode and its damping time for a wider 
range of the flow velocity and for different values of the thread density. We see that for 
a specific value of the flow velocity the parallel fast mode phase speed coincides with that 
of the external leaky slow modes, that is in the range K (ctc) , 3^ ' f^'^ such 

flow velocity, the parallel fast kink wave couples with th e external sl o w mo des by means of 
a "weak" coupling, according to the nomenclature from ISoler et al.l (120081 ). This coupling 
causes a minimum in the damping time as is clearly seen in Figure ^p. On the other hand, 
a similar argument can be adopted to explain the increase of the anti-parallel fast mode 
damping time with f/p, since the phase speed of the anti-parallel wave moves away to the 



-3? Ae ,-3ft(£Tc 



as the 



lich the coupling takes place 



phase speed region of (anti-parallel) external slow modes, i.e. 
flow speed increases. For parallel waves, the flow velocity for w^ 
depends on the thread density, in such a way that the smaller the density, the smaller the 
flow velocity. Hence, the minimum of the damping time moves to smaller values of the flow 
velocity for a small thread density. Nevertheless, this minimum takes place at a flow velocity 
of about 40 km s~^ for a representative thread density of Pp = 5 x 10"^^ km s"^, which 
is a velocity larger than those observed by at least a factor of two. Moreover, even the 
smallest damping time is several orders of magnitude larger than the lifetimes of prominence 
threads, meaning that the effect of the flow is not enough to obtain a reasonable and realistic 
attenuation for kink modes. 
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Fig. 8. — a) Phase speed and b) damping time versus the flow velocity for the parallel fast 
kink mode with kza — 10~^. Different linestyles represent different values of the prominence 
density: pp = 5 x 10"-^^ kg m~^ (sohd line), pp = 4 x 10"^^ kg m~^ (dotted hne), pp = 
3 X 10~^^ kg m~^ (dashed line), and Pp = 6 x lO"^"*^ kg m^'^ (dot-dashed line). The shaded 
zone in panel a) shows the region of phase speeds corresponding to external leaky slow modes. 
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3.2.3. Non-adiabatic sound speed 



Next, we study the influence of tlie steady mass flow on tlie value of the internal non- 



adiabatic sound speed, Ap (eq. [8]). Since 3? y^pj corresponds to the value of the flow 
velocity for which anti-parallel slow waves become backward waves, it is interesting to assess 
the behavior of 3? ( Ap ) as a function of the wavenumber and the flow velocity. We perform 



a similar study as in lCarbonell et al.l (120081) and compute the real part of Ap as a function of 



kz for different values of the flow velocity (Fig. [9]). The behavior of 3? ( Ap j with kz is the one 



explained by lCarbonell et al.l (120081 ) . and presents three different plateaus. For large k^ the 
dominant mechanism is prominence thermal conduction, and the non-adiabatic sound speed 
coincides with the isothermal value. For intermediate kz the non-adiabatic sound speed does 
not depend of non-ideal effects, and its value corresponds to the adiabatic sound speed. 
For small kz, including the observed region of wavelengths, the non-adiabatic sound speed 
becomes slightly smaller than the isothermal one, and prominence radiation is the governing 
mechanism. On the other hand, one can see that the effect of the flow is to widen the range of 
kz in which the non-adiabatic sound speed matches the adiabatic value: the larger the flow, 
the wider the range. For large flow velocities, the transition between the intermediate- /c^ 
(adiabatic) and the small-fc^ plateaus takes place within the observed wavelengths region. 
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Fig. 9. — Real part of the internal non-adiabatic sound speed (in units of the internal 
adiabatic sound speed) versus the longitudinal wavenumber. a) Result for different flow 
velocities: Up — Oknis"^ (solid line), Up — 5kms~^ (dotted line), Up — lOkms"^ (dashed 
line). Up — 20kms~^ (dot-dashed hne), and Up — 30kms~^ (three dot-dashed hne). b) 
Result for Up — lOkms"^ considering all non-adiabatic mechanisms (dashed line), neglecting 
prominence conduction (asterisks), and neglecting prominence radiation (diamonds). The 
shaded zone in both panels indicate the range of observed wavelengths. 
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4. Conclusions 

In the present work, we have studied the combined effect of both non-adiabatic effects 
and a steady mass ffow on the damping of oscillations supported by an individual, cylindrical, 
and homogeneous prominence thread. Our main conclusions are summarized next. 

1. In the absence of flow, slow modes are efficiently damped by non-adiabatic effects, 
while fast kink waves are in practice non-attenuated since their damping times are 
much larger than typical lifetimes of filament threads. 

2. The damping by non-adiabatic mechanisms of transverse, kink oscillations is much less 
efficient in the present, cylindrical case than in the slab geometry. 

3. The thermal wave is a non-propagating solution in the absence of flow and its attenu- 
ation by radiative losses is extremely quick. 

4. The presence of flow breaks the symmetry between waves propagating parallel or anti- 
parallel to the flow. For a flow velocity larger than the internal non-adiabatic sound 
speed anti-parallel slow waves become backward waves. 

5. In the presence of flow, the thermal mode behaves as a wave that propagates parallel 
to the flow and its motions are mainly polarized along the longitudinal direction. 
Nevertheless, this oscillatory behavior cannot be hkely observed in practice due to its 
quick damping. 

6. The damping time of both slow and thermal waves is not affected by the flow. On 
the contrary, for reahstic values of the flow velocity, the larger the flow, the larger the 
attenuation of parallel fast kink waves, whereas the contrary occurs for anti-parallel 
fast kink solutions. Nevertheless, this effect is not enough to obtain realistic damping 
times in the case of fast kink modes. 

In agreement with previous studies, the consideration of non-adiabatic mechanisms pro- 
vides with damping times that are compatible with observations in the case of slow modes, 
which can be related to long-period oscillations. The main effect of the flow on these solu- 
tions is that only parallel propagation to the flow is allowed for strong enough flows. The 
negligible attenuation of fast kink modes in the absence of flow is slightly improved in the 
case of parallel waves when flow is present, the damping time being diminished by an order 
of magnitude for realistic flow velocities. However, the damping time is still several orders of 
magnitude larger than the lifetimes of filament threads and, therefore, neither non-adiabatic 
mechanisms nor mass flows provide with reasonable fast mode damping times applicable to 
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prominences. For this reason, it is likely that another damping mechanism is responsible for 

a more efficient attenuation of transverse thread motions, resonant absorption being a good 
candidate which should be investigated. On the other hand, the investigation of the effect 
of flow, and in particular counter-streaming flows, on the damping of collective transverse 
oscillations of cylindrical multithread models is interesting in the light of the present results 
and could be the subject of a future research. 

The authors acknowledge the financial support received from the Spanish Ministerio 
de Ciencia y Tecnologfa and the Conselleria d'Economia, Hisenda i Innovacio under Grants 
No. AYA2006-07637 and PCTIB-2005GC3-03, respectively R. Soler thanks the Conselleria 
d'Economia, Hisenda i Innovacio for a fellowship. 



Expressions for the perturbations as functions of the divergence of the perturbed veloc- 
ity, A = V ■ vi, and its derivative are given next. 



A. Expressions for the perturbations 
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02 _ p A 2 



= zBo^ -3 ^A, (A9) 



"0 

^0 ~ ^z-^O 



Pmi = ipoVA 7^ ^, (AlO) 



PTi =Pi+ Pmi = ipo F^-^ A. (All) 



B. Region of non-existence of the thermal mode 



Following a similar argument as in Ivan der Linden fc Goossend (Il99ll ) , let us consider 



that the thermal mode frequency in the absence of flow is given by a; = is, with s a real 
quantity. So, from equation (fT6|) one obtains 

^ ^ --{s' + klvl) ^. (Bl) 



rrin 



{vl + Ki){s' + klc%^) ' ^ 

Quantities A and B are the following third-order polynomials in s, 

A = - M2S^ + ky,s - klclN2, (B2) 
B = - Afss^ + ky^s - ky^ATi, (B3) 



with 



(7 - 1) fT, 



= ^ ^-^ {^K\\kl + ujT-ujA , (B4) 

Af2 = - 1) I^^^K^ + (B5) 

ATs = ^^4i4i. (B6) 

vi + ci 

The condition > implies that sign (^p) ^ sign(i3p). The solutions of ^ = are a 
pair of complex conjugate roots and a real root, while the same stands for the roots of 
B = Q. Then, the condition > is only verifled in the region between the real roots of 
Ap = and Bp = 0, namely s_Ap and s^p, respectively, which are very close to each other. 
On the other hand, the external evanescent requirement (m^ < 0) is verified outside the 
region between the real solutions of = and Be = 0, namely and s^^ respectively. By 
computing these real roots (Fig. [TOl) . one obtains that both regions do not overlap except for 
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5.04 X 10~^ < kztt < 5.10 X 10~^. Thus, outside this extremely narrow overlapping region, the 
thermal mode exists with a damping rate in the range sq^ < s < sj^^ where both conditions 
mi > and ml < are satisfied. 
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Fig. 10. — Damping rate per wavenumber versus the wavenumbcr, both in dimcnsionless 
units. The shaded zone points the region where the thermal mode does not exist as an 
evanescent-hke solution, i.e. the overlap of the regions where > (between solid lines) 
and ml > (between dashed lines) are both satisfied. 
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